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Abstract 

In this paper we study the semileptonic decays of the B c meson in the Light-Cone Sum Rule 
(LCSR) approach. The result for each channel depends on the corresponding distribution amplitude 
of the final meson. For the case of B c decaying into a pseudoscalar meson, to twist-3 accuracy 
only the leading twist distribution amplitude (DA) is involved if we start from a chiral current. 
If we choose a suitable chiral current in the vector meson case, the main twist-3 contributions 
are also eliminated and we can consider the leading twist contribution only. The leading twist 
distribution amplitudes of the charmonium and other heavy mesons are given by a model approach 
in the reasonable way. Employing this charmonium distribution amplitude we find the cross section 
cr(e + e~ — > J ftp + rjc) — 22.8 fb which is consistent with Belle and BaBar's data. Based on this 
model, we calculate the form factors for various B c decay modes in the corresponding regions. 
Extrapolating the form factors to the whole kinetic regions, we get the decay widths and branching 
ratios for various B c decay modes including their r modes when they are kinematically accessible. 
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I. INTRODUCTION 

The B c meson has been observed by the CDF and DO groups in different channels 

flag. 

The semileptonic decays of B c was studied in Ref. [4| using the BSW(Bauer, Stech, Wirbel) 
model [5I and the IGSW(Isgur, Grinstein, Scora, Wise) model {(J, and in the frame work 
of the Bethe-Salpeter equation in Ref. [z| and in the relativistic constituent quark model in 
Ref. jg|. Alongside the small differences in the partial decay widths in these models, the first 
estimates made on the basis of the three points (3P) QCD sum rules (SR) [9| are significant 
smaller. The reason was supposed to be the valuable role of Coulomb corrections, which 
implied the summation of a s /v corrections significant in the B c [k]]. It is suggested that the 
discrepancy observed between the QCD sum rules and the quark models can be eliminated 
by including these higher QCD corrections. 

However, the 3PSR inherits some problems when describing heavy-to-light transitions, 
the main one being that some of the form factors have a nasty behavior in the heavy 
quark limit [111 ]. The reason is, when almost the whole momentum is carried by one of the 
constituents, the distribution amplitude of the final meson can be described by the short- 
distance expansion. Moreover, the calculation for the form factors is valid only at the point 
q 2 = 0, and a pole approximation has to be employed to study the semileptonic decays. 
These limit the applicability of QCD sum rules based on the short- distance expansion of 
a three-point correlation function to heavy-to-light transitions and calls for an expansion 
around the light-cone, as realized in the light-cone sum rule approach. In this paper we 
will try to study the semileptonic decays of the B c meson in this approach and compare the 
results with the traditional sum rule approach. 

The semileptonic decays of the B c meson involve the transition B c — > r] c , J/ip, D, D*, B, 
B* , B s , B*. For the case of B c decaying into a pseudoscalar meson, to twist-3 accuracy only 
the leading twist distribution amplitude (DA) is involved if we start from a chiral current. If 
we choose a suitable chiral current in the vector meson case, the main twist-3 contributions 
are also eliminated and we can consider the leading twist contribution only. The result 
depends on the corresponding distribution amplitude (DA) of the final meson. We have to 
construct realistic models for describing the heavy quarkonium and other heavy mesons. In 
particular, the behavior of rj c and J/ip DA's is an interesting subject by the Belle result 
for the cross section a(e + e~ — > J/ip + rj c ). Hence we pay more attention to discussing the 
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heavy quarkonium DA. We calculate the cross section cr(e + e~ — > J/ip + r] c ) by employing 
our charmonium distribution amplitude and the result is consistent with the experiment 
data. Based on a phenomenological model for the leading twist DA, we calculate the form 
factors for various B c decay modes in the corresponding regions. Then we extrapolate the 
form factors to the whole kinetic regions, and get the decay widths and branching ratios for 
various B c decay modes including their r modes when they are kinematically accessible. 

This paper is organized as follows. In the following section we derive the LCSRs for the 
form factors for various B c decay modes. A discussion of the DA models for charmonium and 
other heavy mesons is given in section III. In section IV the cross section <r(e + e~ — > J/ip+r] c ) 
is calculated by using our charmonium distribution amplitude. Section V is devoted to the 
numerical result for the semiletonic B c decays and comparison with other approaches. The 
last section is reserved for summary. 

II. LCSRS FOR THE B c SEMILEPTONIC FORM FACTORS 

According to the definition, the weak transition matrix element B c — >P{V) can be 
parametrized in term of the form factors in the following way: 

<P(p 2 )\qi,Q\B c ( Pl )> = f + (q 2 )( Pl + P 2), + f-(q 2 )q„ (i) 



< V(p 2 )\q^(l - lh )Q\B c { Pl ) >= -ie*Jm Bc + m v )A 1 (q 2 ) +i( Pl + p 2 )^e*q) 



A + (q 2 ) 
m Bc + my 



+ ^( e g) ~ +m + q ^ + P^ \ + m > ( 2 ) 

m Bc + my m Bc + my 

where q = p\ — p 2 is the momentum transfer, e* is the polarization vector of the vector 
meson. 

For B c — > PW we follow Ref. [12J and consider the correlator ILj(p, q) with the chiral 
current, 



n M (p,g) = % yrf 4 xe^<P(j9)|T{g(x) 7At (l + 75)gi(x),gi(0)^(l + 75)Q2(0)}|0> 

= U + (q 2 ,(p + q) 2 )(2p + q) fl + U_(q 2 ,(p + q) 2 )q ll . (3) 

A standard procedure, concentrating on IT + (g 2 , (p + g) 2 ), results in the following LCSR for 

U(q 2 )-- 

f + (q 2 ) = Mrni + rn 2 )f P emk/M2 ,i m 2 - u{g 2 _- um 2 P ) ^ ^ ^ 

m B J Bc Ja p u uM 2 

(4) 



with u = 1 — u and 

A P = [ v /(^- g 2- m 2,)2 + 4m 2 >(m 2_ 9 2) _ ( fl P _ g 2 _ m | )]/( 2m | ), ( 5 ) 

where mi is the mass of the decay quark Qx, m 2 the mass of the spectator quark Q 2 , and 
So and M 2 denote the corresponding threshold value and the Borel parameter respectively. 
In deriving Eq. (j4j) the following definition of the leading twist distribution amplitude (DA) 
<p(u) of the pseudoscalar meson has been used: 

<P{p)\Tq(x)^ 5 Q(0)\0> = -iVp.jp [ due iupx tp(u) + higher twist terms, (6) 

J 



with u being the momentum fraction carried by q. It has been pointed out in Ref. [12j that 
all the twist-3 contributions have been eliminated so those DA's entering the higher twist 
terms in Eq. (J3J) are of at least twist 4. By repeating the procedure for II_(g 2 , (p + q) 2 ) we 
find a simple relation between f+(q 2 ) and f-(q 2 ) up to this accuracy: 

f-(Q 2 ) = -hid 2 ) (7) 
For B c — ► Vlv we choose the following correlator as our starting point: 



n„(p,g) = -i J rf 4 xe^<V(p)|r{g(x)7^(l-75)Qi(x),gi(0)(l+75)Q 2 (0)}|0> 

= Y% - T + (e*q)(2p + q)q, - T-(e*q)q, + iT v e^ fh e* a (fp\ (8) 

Also we take the standard definition of the twist-2 and twist-3 distribution amplitudes of 
the vector meson (see, e.g., Ref. [131]). and neglect higher twist DA's which is supposed to 



be less important in comparison with those written below: 

< V(p)\qp(x)Q a (0)\0 > =y\ U e i ^{f v m v [e*g ( l\u)+p^( ( f> l] (u)- g^(u))} 



- ifW^P'Mu) + ~±Uv ~ fv— ^K a ^e*»p a xVgf>(u)} a pl9) 

In Eq. (Q u is also the momentum fraction of q, and m g (Mg) is the mass of q(Q). 

Similarly one can obtain the following sum rules for Ai(q 2 ), A±(q 2 ) and V(q 2 ) in Eq.(2): 

Mq 2 ) = fUrni + m 2 ) %JM , 
f Bc m Bc {m Bc + m v ) 



1 du 
— exp 

A v U 



m\ — (1 — u)(q 2 — urriy) 



uM 2 



m 2 -q 2 + u 2 m 2 v 

<P±.\u), (10) 

u 
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with 



A 



/ 2n /v(mi+m 2 )K+m F ) m 2 /Jvf2 

-w J = -r~~2 e c 



fB c mj 



1 (iit 
— exp 

A v U 

2 



m\ — (1 — u)(q 2 



urn 



v) 



uM 2 



A-(Q 2 ) = -A + {q 2 ).. 
V(q 2 ) = A + {q 2 



Ay = [ v /(^-g2_ m 2 /) 2 +4m 2 /(m 2_ g2) _ {g V _ g 2 _ ^y^m'y), 

and also mi the mass of the decay quark Q\, m 2 the mass of the spectator quark Q 2 . 



(11) 
(12) 
(13) 

(14) 



III. THE DISTRIBUTION AMPLITUDES OF THE CHARMONIUM AND 
OTHER HEAVY MESONS 

The leading twist distribution amplitude for the heavy quarkonium, as defined in the 
previous section, can be related to the light-cone wave function i/j m (x, kj_) as: 



where J'm is the decay constant. In the non-relativistic case, the distribution amplitude 
<Pm{x) goes to the 5-like function and the peak is at the point x = 1/2. For heavy quarko- 
nium, r] c , the DA should be wider than the 5-like function since the c quark is not heavy 
enough. Of course, it goes to 5-function as the heavy quark mass m* — > oo. 



For the massive quark- ant iquark system, Ref. [bjj provides a good solution ip, 



Aexp^—b 2 ^ 2 ) of the bound state by solving the Bethe-Salpeter equation with the the har- 
monic oscillator potential in the instantaneous approximation. Then one can apply Brodsky- 
Huang-Lepage (BHL) prescription [151 ] : 



i>c.M(q 



J2\ 



LC 



'k 2 ± + m* 2 



^ x(l — x) 

and get the momentum space LC wave function: 



— M 



ip M (x,k ± ) = A M exp 



-b 



k 2 , + m 



*2 



M 



(16) 



(17) 



x(l — x) 

where rriq is the heavy quark mass and M is the mass of the quarkonium. Furthermore, the 
spin structure of the light-cone wave function should be connected with that of the instant- 
form wave function by considering the Wigner-Melosh rotation. As a result, the full form of 



the light-cone wave function should be 

ip f M (x, k±) = xm(x, k±)ip M (x, kj.) (18) 

with the Melosh factor 

XM(x,k ± )= ™ Q (19) 

V k i + m h 

After integrating out k^, the leading-twist distribution amplitude of the heavy quarkonium 
becomes 

8vr a / } M b M ^ x (i _ x ) 

where Erf(x) = ~ f$ exp(—t 2 )dt. As rriq — > oo, <^m(x) goes to the 5-like function certainly. 
This model of the rj c distribution amplitude has been used to study the large-Q 2 behavior of 



rjc-j and ^-7 transition form factors in Ref. 161 ] . The parameters Am and b 2 M in Eq. (fTTj) can 



be determined by two constraints on them completely. One constraint is from the leptonic 
decay constant f M 

(21) 

and another one from the probability of finding the \QQ > state in the heavy quarkonium, 

1 dx I ^l?M^k ± )| 2 = P M . (22) 

— 1 for heavy quarkonium. Inputting the constituent mass m* ~ 1.5 GeV, and the 
decay constant f Vc ~ 0.40 GeV 1 , We get the corresponding parameters for r\ c : 

A Vc = 128.1 GeV, b Vc = 0.427 GeV" 1 (23) 



Then the behavior of the leading twist ru DA can be given and the 



comparison with the 



model from the QCD sum rule analysis 17] and the model in Ref. 13| is plotted in Fig. [IJ 
The moments of these models are given in Tab. [B All the wave functions and corresponding 
moments are defined at the soft scale /x* ~ lGeV. However, the appropriate scale fi for the 
wave functions entering the LC sum rules will be fi ~ rrib for 6-quark decays and fi ~ m c for 
c-quark decays with m h and m c the one loop pole masses. Since /1 is not far from this 



1 The value of /m is taken from the leptonic decay of J/tp: F(J/ip — > e+e~) = (IQira 2 /27)(\fj/^\ 2 /Mj/^), 
fj/^p ~ 0.41GeV. The one loop corrections (~ a s /ir) to the ratio T(rj c — > 2 r y)/T(J/ijj — ► e + e~) indicate 
that f Vc is slightly smaller than fj/^ and we take f rjc ~ fji^ — 0.40GeV on average. 
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TABLE I: The moments of our model for the r\ c distribution amplitude, compared with that in 
Ref. Q and Ref. Q- 



<e n > 


This work 


[13] 


[12 


n=2 


0.21 


0.13 


0.07 


n=4 


0.053 


0.040 


0.012 


n=6 


0.018 


0.018 


0.003 



scale dependence can be neglected in our calculations for simplicity. From Tab. [I] it can be 
found that the moments of the model (20) is similar to that in Ref. [kJ, but much larger 
than that in Ref. 171 ] . Obviously the Melosh factor Xm{x, kj_) — > 1 in the heavy quark limit 
m*n — ■> oo. If we neglect this factor and integrate from Eq. (fT7"|) . we get the corresponding 
distribution amplitude which has a much simple form: 

VSA M 



-x(l — x) exp 



x(l — x) 



(24) 



Sn 2 f M bl 

Actually this is just the wave function proposed in Ref.[l7j based on the QCD sum rule 
analysis. 
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FIG. 1: The model for the leading twist distribution amplitude for ry c (in solid line), in comparison 
with the one in Ref. Q (dashed line) and Ref. [17( (dash-dotted line). 
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For the vector charmonium, J/ip, it is expected that the behavior of the transverse 
distribution amplitude is the same as that of the longitudinal DA since there is no light 
quark in the charmonium system, ie.: 

0||Or) = <MaO = <(x) (25) 

which is confirmed by the moment calculation in the QCD sum rules fis| . 

For the D, B and B s meson, which are composed by one heavy (Qi) and one light quark 
(Q2), According to the BHL prescription one takes the following connection: 



"0c.m(<? 



7?2\ 



LC 



'k\ + mf + kl + mf 



Mi 



(26) 



x 1 — x 

with m\{m^) the constituent quark mass of Qi(Q2), x the momentum fraction carried by 

Qi. Also the Melosh factor should be modified as 

. , , (1 — x)rri\ + xm% 

X p{x, kj.) = v 2 

Jkj_ + ((1 — x)m* + xml) 2 
From which we get the light-cone wave function for pseudoscalar meson 



(27) 



ipUx,kj 



A p xp{x, k ± )exp 



2 (k\ + mf t kl + mf 

-On 



+ 



X 



1 — X 



(28) 



and the corresponding distribution amplitude 2 
V6A P y 



ipp(x) 



%^/ 2 fpb P 



x(l -x)[l -Erf( 



b P y m _ r l2 (xmf + (1 - x)mf - y 2 ). 



x(l 



=)]exp[-6 P - 



x(l — x) 



(29) 

where y = xm* 2 + (1 — x)m\. Similarly, there are two constraints Eq. (f2~Tj) and Eq. (|22il to 
determine the unknown parameters. We take Pd — 0.8, Pp — Pb s — 1.0 as suggested 



in Ref. 2l|]. Inputting the decay constants 
results reported by the CLEO Collaboration 
27 



We use the least-squares fit values o 



22] and lattice simulations 
281 ] . see Tab. [TT]) and the constituent quark masses m* = 0.35 GeV, m* = 



23 



24 



25 



the 



26j, 



0.5 GeV, ml 



1.5 GeV, m£ = 4.7 GeV, we get the parameters: 

A D = 116 GeV b D = 0.592 GeV" 1 , 

A B = 1.07 x 10 4 GeV b B = 0.496 GeV" 1 , 

Ap s = 2.65 x 10 4 GeV b Bs = 0.473 GeV" 1 . 



(30) 



2 This model has been used in Refs. 



19 



20l | for the D meson distribution amplitude with the different 



parameters. There was a misprint of the factor \/2 with the decay constant in Ref. [l9| 
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TABLE II: Leptonic decay constants (MeV) used in the least-squares fit for our model parameters. 





This work 


other 


Sd 


223 


222.6 ± 16.7±|;| CLEO[22j 






201 ± 3 ± 17 MILC LAT [23] 






235 ± 8 ± 14 LAT[24J 






210 ± W±H UKQCD LAT [25] 






211 ± 14lf 2 LAT[26j 


Sb 


190 


216 ± 9 ± 19 ± 4 ± 6 HPQCD LAT[27J 






177 ± 17tl UKQCD LAT[25j 






179 ± 18+^ 4 LAT [26] 


fB 3 


220 


259 ± 32 HPQCD LAT[27J 






204 ± 16±l 6 LAT [26] 






260 ± 7 ± 26 ± 8 ± 5 LAT[28] 






204 ± 12±|f UKQCD LAT[25] 



The distribution amplitudes of these heavy-light mesons are plotted in Fig. [2 The distri- 
bution amplitudes of the corresponding vector mesons are treated in the same way as J/ip. 



IV. THE CROSS SECTION a{e + e~ -> J/if> + rj c ) 

Following Ref. the cross section cr(e + e _ ^ J/ip + rj c ) can be calculated by using 
the distribution amplitudes (1201) and (1251) . We neglect the complicated but slow logarithmic 
evolution of wave function forms, and account only for the overall renormalization factors of 



the local tensor and pseudoscalar currents and for running of the quark mass, as in Ref. 13]. 
One obtains 

a(e + e- -> J/ip + rj c ) ~ 22.8 fb (31) 
which is consistent with Belle and BaBar's measurements 
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301 ] of this cross section: 



a(e + e- -> J/ip + r] c ) = 25.6 ± 2.8 ± 3.4 fb (Belle), 
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FIG. 2: The leading twist distribution amplitudes for heavy-light pseudoscalar mesons. 

a(e + e~ -> J/tp + r} c ) = 17.6 ± 2.8±|;f fb (BaBar). (32) 



The value given by Eq. (}3~Tj) is the same order of the numerical result in Ref. [13j and much 
larger than the standard non-perturbative QCD (NRQCD) calculation. The reason is that 
the DA behavior of the charmonium in our paper and Ref. [13j is much wider than 5-like 



function due to the relativistic effect. Also our result confirms the observation by Ref. [31 ]. 
It may be expected that the large disagreement between the experimental data and the 
standard NRQCD calculation can be resolved by combining the light-cone wave function 



with relativistic effect and radiative corrections 



32]. 



V. NUMERICAL RESULT FOR SEMILEPTONIC B c DECAYS 

For the decay constant of the B c meson, we recalculate it in the two-point sum rules 
using the following correlator 

K(q 2 ) = i J d*xe iqx < 0|c(z)(l-7 5 )&(aO,6(0)(l + 7 5 )c(0)|0 >, (33) 

for consistency. The calculation is performed to leading order in QCD, since the QCD 
radiative corrections to the sum rule for the form factors are not taken into account. We 
also neglect the higher power correction corresponding to the gluon condensates. The value 
of the threshold parameter s is determined by requiring the experimental value of B c be 
obtained in the reduced sum rule after taking the derivative of the logarithm of the SR with 
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respect to 1/M 2 . The quark mass parameters entering our formulas are the one-loop pole 
masses for which we use m& = 4.7 GeV and m c = 1.3 GeV (cf. Tab. 3 and Tab. 4 in the 



review 



33| and references therein). To get the experimental value = 6.286 GeV 



find Sq should be s n ^ 42.0 GeV 2 , which is smaller than the threshold value input in the 



ordinary sum rule 34j. This will ensure in some sense that the scalar resonances will make 



less contribution in our sum rule. The corresponding value of /s c is fs c = 0.189 GeV, which 



is smaller than that in Ref. 34j since we do not include the a s corrections. The same set 
of parameters will be used in the LCSRs for the form factor in order to reduce the quark 
mass dependance. Take the derivative of the logarithm of the LCSR for the form factors 
with respect to 1/M 2 , we get a sum rule for the mass of the B c meson . Requiring this sum 
rule to be consistent with the experiment value at q 2 = 0, we can determine M 2 for each 
LCSR. This results in M 2 (B C -> Vc ) = 25.8 GeV 2 ,M 2 (5 c = 11.6 GeV 2 ,M 2 (5 c -> 

B) = 112 GeV 2 ,M 2 (5 c -> B s ) = 111 GeV 2 . It seems that the Borel parameters for 
B c — ► B(B S ) are somewhat large. However, the LCSRs are quite stable in the large region 
of the Borel parameter 50 GeV 2 < M 2 < 150 GeV 2 actually and we just use the above 
value for explicit calculation. For the vector meson we simply use the same M 2 as the 
corresponding pseudoscalar meson just as we do for the DA's. Also we make the assumption 
that f$ = f v = f P . 

With all the parameters chosen, we can proceed to calculate all the form factors involved. 
The results of the form factors at q 2 = are given in Tab. IHIl in comparison with those from 
other approaches. Notice that in our calculation we always have: 

f + (q 2 ) > 0,/_(g 2 )<0,A 1 (g 2 )>0 

A + (q 2 ) > 0,A4q 2 )<0,V(q 2 )>0. (34) 

In the 3PSR approach the same relations can be obtained, but only in the case of nonrela- 
tivistic description for both initial and final meson states, eg., B c — > J/ip(r] c ). In these decay 
modes the QM results show the same signature pattern, as can be seen in Tab. IHII 

Our calculations for the form factors only valid in limited regions where the operator 
product expansion (OPE) goes effectively. For 6-quark decays, the LCSR is supposed to be 
valid in < q 2 < m 2 — 2mfeAQ C D — 15 GeV and for c-quark decays < q 2 < m 2 — 2m c AQ C D — 
0.4 GeV. It turns out that the calculated form factors for can be fitted excellently by the 
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TABLE III: The values of the form factors at q 2 = in comparison with the estimates in the three 
joints sum rule (3PSR) (with the Coloumb corrections included) [l(| and in the quark model (QM) 

A. 



Mode 




/+(0) /-(0) 


Ai(0) A+(0) A_(0) V(0) 


B c -» cc[lS] 


This work 
3PSR [10] 
QM [8] 


0.87 -0.87 
0.66 -0.36 
0.76 -0.38 


0.75 1.69 -1.69 1.69 
0.63 0.69 -1.13 1.03 
0.68 0.66 -1.13 0.96 


B c -> Bi* ] 


This work 
3PSR [10] 
QM [8] 


1.02 -1.02 

1.3 -5.8 
-0.61 1.83 


1.01 9.04 -9.04 9.04 
0.69 -2.34 -21.1 12.9 
-0.33 0.40 10.4 3.25 


B c -> BW 


This work 
3PSR [10] 
QM [8] 


0.90 -0.90 
1.27 -7.3 
-0.58 2.14 


0.90 7.9 -7.9 7.9 
0.84 -4.06 -29.0 15.7 
-0.27 0.60 10.8 3.27 




This work 
3PSR [10] 
QM [8] 


0.35 -0.35 
0.32 -0.34 
0.69 -0.64 


0.32 0.57 -0.57 0.57 
0.43 0.51 -0.83 1.66 
0.56 0.64 -1.17 0.98 



parametrization: 

Fl{q2) = l-a l9 2 /< + M9 2 K) 2, (35) 
Extrapolate the calculated form factors to whole kinetic region using this parametrization, 

we can proceed to calculate the branching ratios of the simileptonic decays of B c . The 

results is shown in Tab. IIVI together with those of other approaches, where we have used the 

following CKM-matrix elements: 

V cb = 0.0413, V ub = 0.0037, 

V cs = 0.974, V cd = 0.224. (36) 

For the 6-quark decay modes in the B c meson, our results for the branching ratios are much 
larger than the corresponding results in the 3PSR approach. In these decays the kinetic 
region is rather large, so the branching ratios depend slightly on the absolute value of the 
form factors at q 2 = 0. In the LCSR approach, the form factors always increase much faster 
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TABLE IV: Branching ratios (in %) of simileptonic B c decays into ground state charmonium states, 
and into ground charm and bottom meson states. For the lifetime of the B c we take t(B c ) = 0.45ps. 



Mode 


This work 3PSR [10] QM [8] 


m 


rj c eu 


1.64 


0.75 


0.98 


0.97 




0.49 


0.23 


0.27 


— 


J/i/jev 


2.37 


1.9 


2.30 


2.30 


J /ipTV 


0.65 


0.48 


0.59 


— 


Deu 


0.020 


0.004 


0.018 


0.006 


Dtv 


0.015 


0.002 


0.0094 


— 


D*eu 


0.035 


0.018 


0.034 


0.018 


D*TV 


0.020 


0.008 


0.019 




Bev 


0.21 


0.34 


0.15 


0.16 


B*ev 


0.32 


0.58 


0.16 


0.23 


B s ev 


3.03 


4.03 


2.00 


1.82 


B*ev 


4.63 


5.06 


2.6 


3.01 



than the simple pole approximation required in the 3PSR analysis, which accounts for the 
discrepancy in these decays. For the c-quark decays in the B c meson, where the kinetic 
region is narrow enough, our results are consistent with the 3PSR approach roughly. 

VI. SUMMARY 

The semileptonic decays of the B c meson are studied in the Light-Cone sum rule approach. 
By using suitable chiral currents, we derive simple sum rules for various form factors, which 
depend mainly on the leading twist distribution amplitude of the final meson. A model 
with the harmonic oscillator potential for the light-cone wave function is employed. Special 
attention is payed to the leading DA of the charmonium. It has been found that our model 
is consistent with the QCD sum rule analysis. Also, the moments are found to be similar 

n 

to the model proposed in Ref. [13J. Based on this model, we calculate the form factors 
for various B c decay modes in the corresponding regions. Extrapolating the form factors 



13 



to the whole kinetic regions, we get the decay widths and branching ratios for all the B c 
semileptonic decay modes. For the 6-quark decay modes in the B c meson, where the kinetic 
regions are quite large, our results for the branching ratios are much larger than the 3PSR 
results. For the c-quark decays in the B c meson, they are consistent with each other in 
general. 

It is a crucial point to construct a realistic model for the light cone wave function of the 
charmonium which is not a non-relativistic subject. Based on the solution of the relativistic 
Bethe-Salpeter equation in the heavy quark system, we provide a model in this paper by 
using the BHL prescription and the behavior of the charmonium DA is much wider than 
5-like function which was employed essentially by the approximation of NRQCD. Thus the 
cross section <r(e + e~ — > J/ip + i] c ) can be enhanced considerably and is about 22.8 fb. 
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